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In magnetized plasmas gravitational and electromagnetic waves may interact coherently and ex- 
change energy between themselves and with plasma flows. We derive the wave interaction equations 
for these processes in the case of waves propagating perpendicular or parallel to the plasma back- 
ground magnetic field. In the latter case, the electromagnetic waves are taken to be circularly 
polarized waves of arbitrary amplitude. We allow for a background drift flow of the plasma com- 
ponents which increases the number of possible evolution scenarios. The interaction equations are 
solved analytically and the characteristic time scales for conversion between gravitational and elec- 
tromagnetic waves are found. In particular, it is shown that in the presence of a drift flow there 
are explosive instabilities resulting in the generation of gravitational and electromagnetic waves. 
Conversely, we show that energetic waves can interact to accelerate particles and thereby produce a 
drift flow. The relevance of these results for astrophysical and cosmological plasmas is discussed. 

PACS numbers: 



I. INTRODUCTION 

In empty and flat space-time linear gravitational wave 
(GW) perturbations and electromagnetic waves (EMWs) 
do not interact when propagating in the same direction. 
If, however, there is a background electromagnetic field 
or a background curvature present the two waves may 
couple, interact and exchange energy 0-0. 

Although GWs typically interact weakly with matter, 
it is of interest to consider how GW-EMW interaction 
may be altered by the presence of a magnetized plasma 
- which is a common state of matter in astrophysical 
and cosmological scenarios. The interaction may; lead to 
production or modification of observable EMWs, as con- 
sidered in the Refs. 0-0, energize the plasma by excit- 
ing Alfven and magnetosonic waves 

0-0 

which may be 

of importance in more complex processes such as super- 
novae explosions, as discussed in Ref. 0. Furthermore, 
it may even modify the expected GW signals that are 
currently under attempted detection and will, hopefully, 
provides us with a new window through which the uni- 
verse can be observed. Studying nonlinearly interacting 
waves may also reveal new types of instabilities that can- 
not be found using conventional linear stability theory, 
see e.g. Ref. [To| . 

Wave interactions are most efficient if they are resonant 
(coherent), i.e. if the frequencies satisfy certain match- 
ing conditions and the relative wave phase remains un- 
changed for a long time. We will thus not consider inco- 
herent (non-resonant) wave interaction, in which case en- 
ergy conversion takes place on a much longer time scale. 
Resonant interaction in vacuum or in the presence of a 
background electromagnetic field requires only that the 
frequencies should coincide since both GWs and EMWs 
propagate along null geodesies and therefore have identi- 
cal phase velocity equal to c, the speed of light. Consider- 



ing GW-EMW interaction in a "medium" the occurrence 
of such resonances is more rare because the phase velocity 
for GW is (in the high frequency approximation and to 
a very high accuracy) still equal to c, whereas for EMWs 
this occurs only for particular wave frequencies (or in the 
limit of very high frequencies). For this reason resonant 
GW-EMW interaction has been studied in Ref. |6j for 
the case of multiple EMWs, where two or more EMWs - 
that seperately do not propagate along null geodesies - 
together produce perturbations that are resonant with a 
GW. See also Ref. ^ljj where generation of GWs with 
the same mechanism was considered but for interacting 
sound waves. 

Ref. |12j showed that resonant interaction between a 
GW and a single low- frequency EMW (a "magnetosonic 
wave" , in the terminology of their formalism) can be re- 
alized for "incompressible" relativistic magnetofluids, i.e. 
with sound velocity equal to c, but no calculations re- 
garding the coupling strength where made. In Ref. 
the excitation of magnetosonic waves by GWs was con- 
sidered with particular focus on almost coherent waves. 
In Ref. |13| kinetic theory was used to derive, among 
other things, dispersion relations governing the coupling 
between a GW and an EMW propagating perpendicular 
to a background magnetic field. From the dispersion re- 
lations in Ref. 01 it is clear that the waves are resonant 
if the wave frequencies coincide with the electron plasma 
frequency (defined below) but the strength of the inter- 
action cannot easily be deduced from those results. 

In this paper we study the interaction between GWs of 
small amplitude and EMWs in a magnetized plasma mod- 
elled by multifluid equations. In the case of perpendicu- 
lar propagation with respect to the background magnetic 
field, we take the EMW to be a high-frequency extraor- 
dinary electromagnetic wave (using standard terminol- 
ogy of plasma wave theory), the only wave in this case 
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for which we may have phase velocity equal to c (with- 
out introducing a more complicated background state). 
The main focus of our study is on the case of parallel 
propagation, where we make use of an exact (when ne- 
glecting gravitational effects) EMW solution that can be 
reduced to e.g. high-frequency electromagnetic waves, 
whistler waves, low- frequency Alfven waves or waves in 
electron-positron plasmas. This solution also allows the 
presence of a relative drift flow of the fluids constituting 
the plasma. The inclusion of a drift in the background 
state is important for at least three reasons, when consid- 
ering GW-EMW perturbations. Firstly, it increases the 
number of ways that resonant interaction can occur, sec- 
ondly, it alters the coupling strength between the waves 
and, most importantly, it supplies the system with free 
energy. This later fact means, as we will demonstrate, 
that the background configuration may be unstable, lead- 
ing to simultaneous generation of GWs and EMWs. 

The paper is organized as follows. In section II we 
give an overview of certain principles of resonant wave 
interactions. The basic equations are presented in sec- 
tion III, i.e. the Einstein-Maxwell system together with 
multifluid equations. We demonstrate in section IV that, 
provided the conditions for the high-frequency approxi- 
mation are fulfilled, GWs can be taken to be in the trans- 
verse traceless gauge also in the presence of matter and 
we derive the corresponding evolution equations, describ- 
ing how the GWs are coupled to the matter perturbation. 
In section V we derive evolution equations for the EMWs, 
including the effect of GWs on the EMWs and expressions 
for the wave energy density. It is shown that circularly 
polarized EMWs can have negative wave energy for cer- 
tain background parameter values and simultaneously be 
resonant with a GW. Resonant interaction between the 
GWs and EMWs is studied in section VI. Wave inter- 
action equations are derived and their solutions are ex- 
amined, whereby the characteristic time scales for con- 
version between gravitational and electromagnetic waves 
are found. For the case of negative energy EMWs the 
coupling to GWs give rise to explosive EMW-GW insta- 
bilities or, with different initial conditions, acceleration 
of plasma flow. The results are summarized and further 
discussed in section VII, including the relevance for as- 
trophysical and cosmological plasmas. 



II. OVERVIEW OF RESONANT WAVE 
INTERACTIONS 

The natural wave modes of the gravitational field in 
the high-frequency approximation [l6j and of the elec- 
tromagnetic field in magnetized plasmas |l7j| are well- 
known. Based on these solutions, we choose a set of wave 
perturbations, F n , that we represent with the following 
complex notation: F n = f n + c.c, where /„ = fne 16 ^', 
f n is the complex amplitude, 9 n the (real) wave phase 
and c.c. stands for complex conjugation of the preceding 
term. Letting /i denote the GW (e.g. some linear combi- 



nation of components of the perturbed metric tensor) and 
f 2 , say, the electric field of the EMW. From the govern- 
ing equations, the Einstein-Maxwell-plasma fluid system 
presented in section lTTTI we derive evolution equations for 
these wave perturbations. As can be found from section 
llVl andlVl these equations take the form 

D 1 f 1 =S 1 (f 2 J 1 ) + ... (1) 

D 2 f 2 = S 2 (f u f 2 ) + ... (2) 

where D n is the linear wave propagator and the inter- 
action source term S n (f m , /„) is an algebraic expres- 
sion [14j involving f m (and possibly also /„) such that 
S n = 6„e l< ^™ is resonant with f n , i.e. the relative phase 
O n — <t>n is constant. The dots in Eqs. and @ indicates 
that there are in general also terms that are non-resonant 
with /„. Non-resonant terms are discarded, as they only 
have an effect on a much larger time-scale. The source 
term Si have the structure such that in the absence of 
EMWs, f 2 = 0, Si vanishes and Eq. reduces to the 
free wave equation for GWs, {d? — c^]/i = 0. Conversely, 
if there are no GWs then S 2 vanishes and Eq. (J2J) re- 
duces to the plasma EMW equation that corresponds to 
the assumed perturbation. 

Our immediate purpose is to determine the coupling 
strength and thereby the characteristic time-scales for 
the wave interactions. We do this by studying spa- 
tially uniform monochromatic waves propagating in a 
(locally) static and uniform background (the relevance 
for less idealized situations will be discussed in Sec- 
tion [^nj . The interaction will then result in time de- 
pendent wave amplitudes, governed by the derived evo- 
lution equations. Since the source terms are small - 
either proportional to the gravitational coupling con- 
stant or to the small GW amplitude - we assume that 
Iwnfnl ^ |<9tf„|. This means that, to lowest order, the 
waves will propagate freely according to the dispersion 
relation D n {— iuj n , ik„) = and, to the next level of ac- 
curacy, have a slowly evolving amplitude. Applying this, 
the differential operators in Eqs. and @ can be ap- 
proximated by D ps D(— iu),ik) + (dujD)dt, where the 
tilde notation on d t means that the derivative acts only 
on the amplitude. Higher order derivatives are neglected. 
The evolution equations and J2J thus becomes 

dth = {du, 1 D 1 )- l S 1 {h,h) (3) 
dtf2 = (d ul2 D 2 )- 1 S 2 (f 1 J 2 ) (4) 

and we refer to these as the interaction equations. De- 
termining Si and S 2 and studying the solutions of these 
equations are the main purposes of this paper. 

III. BASIC EQUATIONS 

We take the gravitational and electromagnetic field to 
be governed by the Einstein field equations 

Gab = nT ab (5) 



3 



and the Maxwell equations 



Maxwell and fluid equations can be formulated as Q , 0] 



V a F bc + V b F a 



V a F ab 







(6) 
(7) 



where G ab is the Einstein tensor, k = 8ttG, T ab is 
the energy-momentum tensor, F a i is the electromagnetic 
field tensor, j a is the total four-current density and V 
denotes covariant derivative. We use units where the ve- 
locity of light in vacuum is c = 1 . The metric signature is 
( — h ++) and for the indices we use a, 6, c, ... = 0, 1, 2, 3 
and fc, ... = 1, 2, 3. The matter present in the interac- 
tion region is a magnetized plasma for which we choose 
a multifluid description. This means that we take the 
plasma to consist of a number of interpenetrating charged 
fluids, one for each species of particles that constitutes 
the plasma. The fluids interact only through the electro- 
magnetic and gravitational field, i.e. we neglect effects 
of particle collisions as can be done for most plasmas. 
The appropriate fluid equations arc then, for each fluid 
component (i), the equation of mass conservation 

V a {m {i) n {C) u a {{) ) = (8) 

and the momentum equation 

m (l) n (i) u a {i) \7 a u b {i) = -V h p (l) + q(i)n {i )U a {l) F b (9) 

where is the proper particle number density, u%\ is 
the fluid four-velocity, is the particle mass, p^ is 
the pressure and q^) is the particle electric charge. For 
closure this should be supplemented by some equation of 
state which we do not specify here. The total current 
density is then given by j a = Q(i) n (i) u (i) an< ^ the 

energy-momentum tensor by T ab — T^P^ + T^f m \ where 
the fluid contribution is 



T ab = ^2( m (*) n (i) + P(i)) u {i)aU{t)b + P(i)9ab (10) 

and the electromagnetic field contribution is 



rp(em) _ p CTJ 



habF cd F ( 



(11) 



From now on we will omit the species index, (i), unless 
there can be any confusion. 

It is practical to rewrite the Maxwell and fluid equa- 
tions in a more convenient form. In particular we in- 
tend to introduce an orthonormal frame so that the in- 
terpretation of the quantities becomes more direct. First, 
we note that by introducing an observer four velocity, 
V a , the electromagnetic field can be decomposed, rela- 
tive to V a , into an electric part, E a , and a magnetic 
part, B a , according to F ab = V a E b - V b E a + e abc B c , 
where E a = F ab V b , B a = \e abc F bc , e abc = V d e abcd 
and e abc d is the 4 dimensional volume element with 
e oi23 = v/| det $|. If one chooses an orthonormal frame 
(in which V a = (1, 0)) with contravariant basis {e a }, the 



V ■ E = p + p E 

V • B = p B 
e E- V xB = -j-j £ 
e B + V x E = -] B 

e o(7?i) + V • (7nv) = An 
mn(e + v • V)7v = — 7 _1 Vp 

+ qn(E + v x B) 



mng 



(12) 
(13) 
(14) 
(15) 
(16) 

(17) 



where we have introduced an Euclidean three- vector no- 
tation E = (E 1 ^ 2 ^ 3 ) etc., V = (ei,e 2 ,e 3 ), p = j° 
and coordinate velocity, v a , so that u a = ^fv a with 
7 = (1 — u,;u 4 )~2 . The effect of the gravitational field is 
included in the form of effective charges, currents, gravi- 
tational forces etc. that were introduced above, and they 
are given by 



Pe 
Pb 

3b 



I''. /••' - e ijk r ;;B k 



ijFk 



I '„/;' + e ijk T ^E, 



V 1 



-e* k (T%B k +T™ k B m )] ei 



V A Oj ~ r )o) BJ i" 1 oj J 
+e«* {T° E k +T™E m )]e l 

g = -7 [r l 00 + (r* 0J + T) ) v> + T) k v>v k ] e, 
where T a bc are the Ricci rotation coefficients. 

IV. GRAVITATIONAL WAVE EVOLUTION 

In this section we show that, in the high-frequency 
approximation |16| for GWs propagating in matter, the 
GWs can be taken to be in the transverse trace-less (TT) 
gauge and we give the corresponding evolution equation. 
In flat space-time linearized gravitational waves can al- 
ways be taken to be in the TT gauge. In vacuum space- 
time with a background curvature this is in general only 
possible in the high frequency limit, where the ratio of 
the GW wavelength and the characteristic background 
length scale (radius of curvature) tends to zero. The 
advantages of having GWs in the TT-gauge is that the 
evolution equation becomes more simple, the polariza- 
tion state is more clear and, most important for us here, 
it reduces the amount of algebra when calculating the ef- 
fect on a plasma, e.g. when calculating the gravitational 
source terms in the EMW evolution equation. However, 
the (perturbed) energy-momentum tensor should have 
the same properties as the perturbed Einstein tensor and 
this will in general not be the case if one just assumes 
the TT-gauge. Nevertheless, we demonstrate that the 
TT gauge can also be applied in the presence of matter 
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provided the conditions for the high frequency approxi- 
mation is fulfilled. This result seems to have been over- 
looked in the litterature. 

Let the background gravitational field and the un- 
perturbed plasma fulfill the Einstein field equations, 
G a ^ b = kT^P , and introduce a perturbation so that 

Gab = Gg 5 + SG ab and T ab = + 5T ab . We assume 
the high frequency approximation so that the background 
space-time can be taken to be Minkowski space-time [T^ | 
and put g a b = Vab + h a b, where r\ a b is the Minkowski 
metric and h ab the small metric perturbation that should 
obey 

SGab = KST ab (18) 

We limit ourself to linear gravitational perturbations. 
The perturbed Einstein tensor, linearized in h a b, is given 
by 

SG ab [h] = -\d c d c h ab + d c d {b h a)c - \iia b d c d d h cd (19) 

h ab = h a b — \"Habh and the brackets / ) stands for sym- 
metrization with respect to the enclosed indices. In vac- 
uum, the next step would be to apply the Lorentz gauge 
condition, d b h a b = 0. Instead of this, we split h a b ac- 
cording to 

Kb = hab + fab (20) 

where h^ b is the (maximal) part of h a b that fulfills the 
Lorentz gauge condition and the other part, f a b, con- 
taining terms that cannot be fitted into h^ b (both terms 
should be symmetric) . With this splitting Eq. (|18f) reads 

-\d c d c h L ab + SGabif] = ^T ab (21) 

Next we assume that the metric perturbation (and the 
perturbed energy- momentum tensor) is of the form h a b = 
t)ab(x d )e tkcxC + c.c. with the dispersion relation k a k a — 
0. The dependency of t) a h on x d is assumed weak, i.e. 
\d c t)ab\ < |Ma&|- Computing 6G ab [f] gives to lowest 
order 

SGabif] = -k c k {b fa)c + \vabk c k d J cd (22) 

where terms involving slow derivatives (derivatives on the 
weakly varying amplitude) have been neglected. From 
this expression it follows, for perturbations with four 
wave vector k a = (w,0,0,w), that 5G X2 [f] = SG 2l [J] = 
and SGn[f] = SG 22 [f]. Furthermore, nonzero compo- 
nents of f a b should be of order / ~ kST/uj 2 (which must 
be much smaller than unity in order not to invalidate 
the linearization in h a b) for the field equations i|21|) to be 
fulfilled. 

We now turn to the terms in h^ b . A gauge transforma- 
tion x a ' = x a +£ a with£ Q < 1 and£ a = -i£ a (x c )e ikbX " + 
c.c. alters the metric perturbation by 

h a ' b ' = h ab - Ck b - £ b k a + rj ab i c k c (23) 



where it is understood that the amplitude, £ a (x c ), is 
weakly dependent on x c in the same way as h ab . We can, 
however, neglect slow derivatives on £ a as such terms 
only produce small corrections to the f a b part or, finally 
in the wave evolution equations, terms of second order 
slow derivatives. Consequently, £ a can, as in the vacuum 
case, be chosen as linear combinations of h^ b such that 

h a ' b ' = hf T + f ab (24) 

where h^ T is "in the transverse traceless gauge", i.e. 
hjt-rp = —h^ T = hA and hif T = h?f T = hs- The evo- 
lution equations for hA and hg follows from Eq. I|21l) by 
subtracting the n and the 22 components and by - for 
aestethical reasons only - adding the 12 and 21 compo- 
nents. The result is 

d c d c h A - -k (ST n - ST 22 ) (25) 
d c d c h B = -k(ST 12 +ST 21 ) (26) 

which depends critically on the fact that 8Gi 2 [f] = 
8G 21 [f] = and 5Gu [/] = SG 22 [/] . 

To summarize, for the part of the metric perturba- 
tion that represents GWs propagating in the x 3 -direction, 
hj? T , we have derived evolution equations that describes 
how the waves are coupled to the perturbed energy- 
momentum tensor. Clearly, nonvanishing STn — ST 22 
and STi 2 can potentially act to drive or damp GWs. 
The remaining part of the perturbed metric, given by 
fab, should be included for the sake of consistency. The 
perturbed energy-momentum that is not accounted for 
in Eqs. and l(2f)l) produces some gravitational re- 

sponse, which is given by SG a b[f] and found from the 
remaining components of Eqs. (|21J) (due to Eqs. Ij21|l 
and these are purely algebraic equations). As f a b 
represents the self-gravitation of the energy-momentum 
perturbation it is in general important for the evolu- 
tion of the matter and electromagnetic field, but since 
/ ~ kST/uj 2 it is negligible when we consider waves in 
the high-frequency limit. We conclude that in the high- 
frequency limit and for slowly evolving wave amplitudes 
the GWs can be taken to be in the TT-gauge and the 
metric perturbation h a b = h^ evolves according to Eqs. 

and {2H>- In this limit the energy density carried by 
GWs follows directly from the Landau-Lifshitz pseudo 
energy-momentum tensor [T^ 

£gw = h [(dth A ) 2 + (d t h B ) 2 ] (27) 

When studying the effect of GWs on a magnetized 
plasma it is practical to use tetrad formalism and intro- 
duce an orthonormal frame. A contravariant basis that 
corresponds to GWs in the TT gauge with propagation 
in the x 3 = z direction is given by 

e = t , ei = (1 - \h A )d x - \h B d v , 
e 2 = (1 + \hA)d v - \h B d x , e 3 = d z 



5 



and the gravitationally induced terms in Eqs. I|12|l - 117l) 
are (displaying the nonzero terms only) 



g = -57(1 - y z) [vxd t h A + Vyd t h B )] ei 
- ±7(1 - v z ) [v x d t h B - v y d t h A ] e 2 



3e 



27 [i v l - v 2 y)dth A + 2v x v y d t h B ] e 3 

-§ \{E X - B y )d t h A + (E y + B x )d t h B ] ei 



- i [-(£„ + B x )d t h A + (E x - B„)^/i B ] e 2 
js = -\ [{E y + B x )d t h A - (E x - B y )d t h B ] e 1 

- \ \{E X - B y )d t h A + {E y + B x )d t h B ] e 2 

where we have used d z « — dt- 



(28) 
(29) 
(30) 



V. ELECTROMAGNETIC WAVE EVOLUTION 

In this section we present the evolution equations for 
the EMWs that we are considering. Since the GW four- 
wave vector is assumed to fulfill the dispersion relation 
k a k a = 0, they have phase velocity equal to unity. The 
EMWs that can interact resonantly with the GWs should 
then also have phase velocity equal to unity. For sim- 
plicity, we restrict our study to the case of perpendicu- 
lar and parallel propagation to the background magnetic 
field, respectively. 

The plasma is assumed neutral and uniform, with a 
constant background magnetic field, Bt \. The plasma 



frequency is denoted by u> p (i) '■ 
cyclotron frequency by w c /j) 



(»(») <$)/"*(<) 



) 1/2 and the 

2(i)-B(o)/7(i) m (i), where 
nu\ is the unperturbed proper particle number density. 
We also allow the presence of a background drift flow 
in the direction of the magnetic field. The inclusion of 
the drift flow provides the system with free energy that 
may be transferred to the GWs and EMWs during the 
interaction. 

In magnetized plasmas there are many electromagnetic 
wave modes and their properties are well known for flat 
space-time and for linear (small amplitude) waves (see 
e.g. H3)- in order to obtain the effect of GWs on the 
EMWs we derive their evolution equations from the basic 
equations H12[l - i|17|) . including the effect of the gravita- 
tional field. For later reference, we note that the total 
energy density of an EMW is given by 



£em — E* 



(31) 



where E cx e %[ - kz ~ u ') t is the electric field of the wave, ui 
is the frequency and e is the effective dielectric tensor 
defined such that D = e-E gives the electric displacement 
field. For the cold (i.e. zero pressure) linearized plasma 
fluid equations e is given byj. 



S -iD 
iD S 
OOP 



(32) 



where there is a background magnetic field and drift flow 

\{R + L),D=\ 



v z in the z-direction, S = \ {R + L), D = | (i2 — L) and 



R= 1 



E 

(i) 

E 

(0 

E 

(0 



kv. 



kv. 



(■<)) 



lu 2 {lu - kv z (i) +W e (i)) 



kv. 



LU 2 {UJ - kv z (i) - W c (j)) 



Although this result concerns linear waves in a cold 
plasma it also applies for the large amplitude waves in 
a moderately warm plasma that we will consider below. 
This is because they are circularly polarized and purely 
transverse, making them effectively linear [20j and with 
no density (nor pressure) perturbations. 



A. Perpendicular propagation 

In the case of perpendicular propagation the plasma 
equations are difficult to treat without resorting to per- 
turbation theory. We restrict ourself in this case to lin- 
ear perturbations and to this level of approximation the 
wave perturbations do not produce any effective force, 
g, that is not rapidly oscillating. This implies that the 
waves cannot exchange energy with the drift flow in an 
efficient way. Consequently, we discard the drift flow en- 
tirely and only considers one particular wave mode in an 
electron-ion plasma, namely the extraordinary electro- 
magnetic wave - using standard plasma physics termi- 
nology. This is the only wave that can be resonant with 
a GW in the simplest plasma model, the cold plasma fluid 
model. 

We take the background magnetic field to be B/q) = 
S/ )ei. The extraordinary EMW is a high-frequency 
wave, to 3> uj p i (such that only the electron fluid is per- 
turbed), and has, when propagating in the z-direction, 
the following non-zero linear perturbations: 5n, v = 
v y &2 + v z G3, E = E y e2 + E z e^ and 5T5 — 8B x e\, where 
the density and velocity perturbation referres to the elec- 
tron fluid only. For a cold plasma, the Eqs. (|12|) - (|17|l re- 
duces to the following equations, linearized in the small 
perturbations 



(33) 
(34) 
(35) 
(36) 

(37) 





d z E z 


= qSn 




dtEy- 


d z SB x 


= —qnvy 


" JEy 




d t E z 


= -qnv z 




d t SB x 


- d z E y 


= - 3B X 






OtVy 


= He v 

m 


+ v z B (Q) ) 




d t v z 


= ^{E Z 


-VyB (0) ) 






m 


dtSn H 


- nd z v z 


= 





(38) 
(39) 
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which can be combined into the following evolution equa- 
tion for the wave magnetic field (after one time integra- 
tion) 

D X SB X = S x (40) 

where the wave propagator, Dx, and the source term, 
Sx, are 

D X = df + (u£ + co 2 )df - d?d 2 - uldl + 4 

S X = \B {Q) [a t 4 + (^ +W 2 )a 2 +w 4 

+d z d t (d 2 + u; 2 p + uj 2 )] h A 

and we have made use of Eqs. Q2H |) -Q3U [I . linearized in the 
wave perturbations and introduced the hybrid frequency 
uj\ = LUp + uj 2 . In the absence of GWs Eq. l|30)) reduces 

to Dx5B x = and, in particular, for wave perturbations 
SB X cx e lt v kz ~ ut ) it produces the dispersion relation for 
the extraordinary wave 

D x =uj 4 - (wjj + u 2 p + k 2 )u? + uo 2 h k 2 + uj* = (41) 

It should be pointed out that the extraordinary mode, 
with the given polarization, only couples to Ka and not to 
h b ■ Computing the perturbed energy-momentum tensor 
gives, linear in the perturbations, 

ST n = -8T 22 = -SB X B {0) (42) 
ST 12 = ST 21 = (43) 

The dispersion relation l|41|l implies that the condition 
for having phase velocity ui/k = 1 is lu — ui p . We also 
note that from Eq. JjHJ and Eqs. p5 Jl -l(33 )l it follows, 
after some algebra, that the energy density of a free ex- 
traordinary wave with ui = ui p is given by 

£ x = 2u- 2 J 2 h \5B x \ 2 (44) 



for electron-ion type of plasmas as well as for electron- 
positron type of plasmas. 

We take the background magnetic field to be B( \ = 
_B(o) e 3 and suppose there is a velocity drift, v z , in the 
z-direction. The above mentioned circularly polarized 
EMW, propagating in the z-direction, has the following 
perturbations (of arbitrary amplitude) v = v x e\ + v y e 2 , 
E = E X B\ + E y e 2 and = SB x ei + SB y e 2 , where v and 
<5B are parallel and all perturbations are functions of z 
and t alone. Note that circular polarization implies that 
7 depends only on the wave amplitude and on v z , not on 
the rapidly varying phase. In the following we treat 7 as 
being constant and v z as a constant uniform background 
flow. Small (linear) deviations in 7 and v z can still be 
described (see Section fvTI and rvTT|l . For convenience we 
introduce the variables E± = E x ±iE y , and similarly for 
all other vector variables. In these variables, suitable for 
circularly polarized waves (the plus / minus variables cor- 
responds to the amplitudes of the right/left hand polar- 
ization), those equations of Eqs. I|12l) - (|17l) that governs 
the given perturbations can be rewritten as 



(d t + v z d z )jv± = -^-(E ± T iB v± ± iv z B±) + g± (45) 
m 

d t E± = ±id z B± - Y,qn~fv± - j E ± (46) 

d t B ± = T*d z E ± - JB± (47) 

where the meaning of g± is g± = g x ±ig y and similarly for 
je± and j B ±. Recall that there are two sets of Eqs. 
(|37Jl . one set for each particle species. These equations 
imply the following evolution equation for E± 

D E E ± = S E (48) 
where the wave propagator is 



B. Parallel propagation 

If one neglects general relativistic effects, there exist 
exact solutions of the multifluid equations that represents 
EMWs propagating parallel to a background magnetic 
field (see e.g. 0). These solutions describe circularly 
polarized waves of arbitrary amplitude and arbitrary fre- 
quency and can thus, taking the appropriate limits, be re- 
duced to high-frequency electromagnetic waves, whistler 
waves, low- frequency Alfven waves or waves in electron- 
positron plasmas. The solution can also be extended to 
include a background drift flow [2^ • We derive the evolu- 
tion equation for this wave for a two-component plasma, 
including the effects of GWs. The indices e and i will now 
stand for negatively charged particles (e.g. electrons) 
and positively charged particles (e.g. positive ions), re- 
spectively, but we make no assumptions for the mass ra- 
tio m e /mi. The calculations are therfore equally valid 



D E = GeGiD + uP vi Gji + uj 2 pe Gj e (49) 

U^d 2 -dl 
d{i) = d t + v z{i) d z 
G(i) = d t + v z{i )d z ± iw c (i) 

and the source term is 

S E = -qinG e d t g t ± - q e nGidtg e ± 

± i (uj 2 l v iz G e + tol e v ez Gi - G e G l d z ^j j B ± 

- G e G l d tjE ± (50) 

In the absence of the gravitational source terms the evo- 
lution equation l|48|) reduces to L>eE± = and for this 
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case we in particular note the solution [2l| 

E± = Ee 1 ^'^ 
.k 



B± = ±i-E. 

CO 

v± = iAE± 



A = 



- kv. 



7TO lo{lo — kv z =p LO c ) 



7 = (1 - v±v T - vf) " 



(51) 
(52) 
(53) 
(54) 

(55) 



with the dispersion relation 
= D E = (lo - kv ze =F w ce ) (lo - kv zi =p lo ci ) (lo 2 - fc 2 ) 

- UJ pi (w - kv ze =F Wee) ~ 

- w 2 e (w - fcv zi =F w cl ) (w - fcu ze ) 
that can be written more compact as 



w 2 -fc 2 



(56) 



(57) 



We again emphasize that the solution is valid for arbi- 
trarily large amplitudes. The condition for the plasma 
to be electrically neutral, j° — 0, implies 7 e n e = 7^. 
As 7(i)^(i) is a preserved quantity (by Eq. I|16[l). the 
proper density, n.(j), is not. Accordingly, both and 
i^ c (i) depend on 7^). This makes the dispersion relation 
amplitude dependent. 

For future reference we devote the remainder of this 
section to establish some properties of the free EMWs. 
As we will focus on EMWs with lo — k we note that 
the dispersion relation in this case implies the following 
relation between the background magnetic field and the 
drift velocities 



±(7e + It) — 

m e v ze - v z 



(58) 



where Q ce = j e to ce /uj 



q e B(Q\/m e is the normalized 



(and preserved) "electron" cyclotron frequency and we 
have made use of the neutrality condition "f e n e = 7iTij. 
Note that for any given background plasma parame- 
ters (except the absence of any drifts), there always ex- 
ist a wave frequency fulfilling Eq. I|58(l . The compo- 
nents of the perturbed energy-momentum tensor, for this 
large amplitude circularly polarized EMW, that couples 
to GWs has the values 

6T n = -\(E 2 x +Bl- E 2 y - Bl) + jmnv 2 x (59) 

ST 22 = 1(E 2 + B 2 X - El -BD+J2 imnv 2 y (60) 

ST12 = - (E x E y + B x B y ) + jmnv x v y (61) 

Because the wave is circularly polarized and transversal 
it is effectively of linear nature (the Lorentz force is linear 
in the wave amplitude) and has no density (nor pressure) 
perturbations. Therefore, the wave energy density can 



be calculated by means of Eq. I|31[l . For the given wave 
perturbations, the wave energy density is given by 



£ E = 7 -d^ 2 A ± )\E ± \ 2 
2lo 

w p uj - kv z 

A ± = 1 ~ ~2" 



(62) 



lo 2 lo — kv z =F LO c 



Applying the dispersion relation l|57|l with lo — fc, the 
wave energy density can, after some algebra, be formu- 
lated as 



£k = 



1 



2cj 2 a e a t 



\2io 2 a e a t - Wp e (l - v e + ou) 



-L0 2 pl (l-v, l + a e )]\E ± \ 2 (63) 
where we have introduced the abbreviation 

"(i) = 1 - v {l) T w cW /w (64) 




FIG.l. The dependence of wave energy Ee (normalized 
by \E±\ 2 and resettled) on drift flow v = Vi — —v e for 
null geodesic small amplitude waves in an electron-ion 
plasma (m e /irii ~ 1000,). The different curves, from the 
top and down, corresponds to the plasma frequencies 
n pe = 2, 5, 10, 20, 40, 60, 100, 1000. As the sign of the 
energy is the main focus here, the energy in the different 
cases are rescaled by suitable factors. 



It is important to note here that Ee is, in contrast 
to Eqw and Ex, not positive definite. For some back- 
groundparameters the wave energy density may be neg- 
ative 10]. This is due to the presence of the velocity 
drifts, providing the system with free energy. Negative 
energy waves are wave perturbations such that the to- 
tal energy of the plasma is less than the energy of the 
unperturbed state j^. To show that there are waves 
fulfilling lo — k and Ee < simultaneously, we consider 
an electron-ion plasma (m e /mj ~ 1000) with drift flows 
v e = —Vi. In FIG 1 the wave energy factor Ee/ \E±\ (for 
small amplitude waves) is presented, as a function of ion 
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drift velocity, for a number of different values of the nor- 
malized "electron" plasma frequency fl pe = "f e Lo pe /u> and 
with the cyclotron frequency fulfilling the condition l|58|) . 
As the sign of the energy is the main focus here, the en- 
ergy factor in the different cases are rescaled by suitable 
factors. Negative wave energy occurs for positive ion ve- 
locities |24j . The minimum required ion velocity depends 
inversely on the plasma frequency. For Q pe — 1000, 
negative wave energy requires only Vi > 10~ 4 whereas 
for rip e ;$ 20 it requires highly relativistic velocities, 
Vi > 0.95. 



etc. We will use the notation of "slow derivatives", 
dt, that acts only on the wave amplitude, so that e.g. 
dtfiA = i^G^A + dth-A, and we will neglect slow deriva- 
tives that are of second or higher order. Consistent with 
this level of approximation, we make use of the relations 
between the perturbations valid for free waves, i.e. Eqs 
(|51|I - H53|I . in the computation of the source terms. Cor- 
rections to this approximations produces higher order 
source terms. 

We now treat the two cases of propagation direction 
separately. 



VI. RESONANT WAVE INTERACTION 

The wave evolution equations , (J2U , 1001) and 

where derived with no assumptions on the specific wave 
forms. In this section we assume the wave frequencies 
and the (uniform) background plasma to be in a state 
where wave resonance occurs, since this gives the most 
efficient energy transfer, and we derive the corresponding 
wave interaction equations. We let 

hA,B = h,A,B + c.c. 
SB X = B + c.c. 
E± = E x ± iE y 

with h A , B = t)A,Be* kG " xa , B = <Be lkx <> xa and E± = 
(£e lkE aX and the four-wave vectors k G = (ujg, 0, 0, ka), 
k\ = (u}x,0,0,kx), and 0, 0, ks) are assumed 

to fulfill the dispersion relations k G a k G — 0, Eq. (|41|) 
and Eq. (|57|) . respectively. As before, c.c. stands for the 
complex conjugate of the preceding term. The conditions 
for the EMWs to be resonant with a GW are, in the case 
of perpendicular propagation 

k a G = k a x (65) 

and in the case of parallel propagation 

k G = 2k a E (66) 

Here we have used that the extraordinary EMW con- 
tributes linearly to the perturbed energy momentum ten- 
sor (|42|) and the circular polarized large amplitude EMW 
quadratically in (|59|I - H61I) and thus produces terms pro- 
portional to e l2 i k EZ-UEt) _ Consistently, the source terms 
for the extraordinary EMW is linear (in the GW ampli- 
tude) whereas for the circular polarized large amplitude 
EMW there are terms proportional to E±\ia.b- Since (in 
the high frequency limit) oja/kc = 1, also the EMWs 
should have tox/kx = 1 and iOE/kE = 1 in order to re- 
main in phase with the GW. The interaction is weak in 
the sense that the coupling terms in the evolution equa- 
tions (|2SJ), 601 and igHI are either proportional 
to the small GW amplitude or the gravitational cou- 
pling constant. This leads to a weak time dependence 
in the wave amplitudes f), 58 and (£, i.e. dt\) <C u>g§ 



A. Perpendicular propagation 

From the dispersion relation l|41|) it follows that in 
order for the extraordinary EMW to propagate along 
null geodesies the wave frequencies must be equal to the 
plasma frequency, ui = to pe . Although this is a spe- 
cial case, a GW or an EMW propagating in a slowly 
varying background density may often reach such res- 
onance regions. Assuming resonant waves, we take k G = 
k x = (uj, 0, 0, u). The wave evolution equations and 
(|40ll then imply the following wave interaction equations 
(dropping the tilde notation on the waves and denoting 
\i = Ha and B = 5B) 

B t h = -iC G B (67) 
d t B = -iC x h (68) 

where Cg = ^ _1 kB(o) and Cx = ^ww fc 2 w^B( ). From 
Eqs. (|27|l and l|44|) the total wave energy is £gw + = 
uj 2 K- x \h^ + 2uj- 2 ujI \B\ 2 and it follows from Eq. (JHZI 
and l|68|) that this is a conserved quantity. The inter- 
action equations has the following solution for the wave 
amplitudes 

B = B cos(V^) + iCx^ho sin(^) (69) 
h = h cos(V>*) + iC G i/j~ 1 B sm(ipt) (70) 

where B = B(t = 0) and ho = h(t = 0),ip = VC X C G = 
v /^|w^ 1 w c -B(o) I- The solution shows that the total wave 
energy alternates between the GW and the EMW with 
the frequency ip. In most applications ip is very small 
and it is then meaningful to linearize the trigonometric 
functions in tjj and it is clear that, for time scales smaller 
than ip' 1 , GWs are converted into EMWs as h = CgBqI 
or, alternatively, EMWs are converted into GWs as B — 
Cxhot. It should be noted that in the later case the 
time t n i for the EMW perturbation to reach the nonlinear 
stage, B ~ -6(0)? is typically much smaller than total 
energy conversion time, . 

B. Parallel propagation 

For resonant interaction between EMWs and GWs in 
the parallel case we take k G = 2k% — (2w,0, 0,2w). As 
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the EMWs that we are considering are circularly polarized 
waves, we also introduce variables for circularly polarized 
GWs h± = h A ± ih B - If the GW and the EMW are 
oppositely polarized the coupling vanishes, so from here 
on we assume identical polarization. The evolution Eqs. 
()25J) and Ij26(l combines to 



(-d 2 + d 2 ) h± = ~Y^ Kmnj 2 v 2 ± 



(71) 



and together with the EMW evolution equation l|48|) this 
implies the following interaction equations for the vari- 
ables h = h± and E = E± 



d t h = iC GW E 2 
d t E = iC E E*h 



where 



Cgw 
C e 



V > 2 a 2 



(72) 
(73) 

(74) 



v zi )a e 



+ 7e — A e ^ £ (i 

Qe 



Ote) 



-it) 



pe 



<)]-' 



(75) 



In the derivation we have noted that the effective cur- 
rents, ]e and js, vanishes and used that g± — —57(1 — 
v z )v±dth±. The coupling coefficients can be related to 
the wave energy densities, Eqs. Ij27(l and (|63|) . according 
to 



Cgw — CE GW \h\ 



Ce — C£ , 



\E\ 



where 



C = — 



2u ^ a 2 (i) 



(1 



(76) 
(77) 



(78) 



It is then straightforward to confirm that Eqs. (|72(l and 
(|73[1 implies that the total wave energy density 



K 



1 



-U) 



pi 



(1 



[2u) 2 a e cti 

+ a e )-oj 2 pe (l-v ze + ai )]\E\ 2 (79) 



2uj 2 a e a l 



is conserved. 

The solution to the wave interaction equations (|72() and 
(|73|1 can be given in terms of Jacobi elliptic functions. 
We exploit the fact that the system can be reformulated 
as a three- wave interacting system (see e.g. 0) with 
two identical EMWs and one GW. This is done by taking 
E 2 = E 3 = \E such that Eq. l(75|) can be splitted into 
two (identical, but differently labeled) equations and the 
right hand side of Eq. (|72"|) reads iA-CgwE^Ej,. Next we 



convert the obtained three- wave interaction equations to 
a real system and renormalize the amplitudes such that 
the coupling coefficients become unity. The real system 
is |25j (dropping the tilde notation) 



d t ipi = s 1 ciip 2 i>3 cos(</>) 

dt?p2 = S2C2l/'lV ; 3 COs(</>) 

d t 4>3 = saCaV'iV^ cos(</>) 



(80) 
(81) 
(82) 



where we have taken ih — tpie 1 ^ 1 , E n — ip n e l( '' n , 
-ACgw = sici, C E = s n c n (for n= 2,3) with tpi = \h\, 
4>n = \E n \, ci = I Cow - 1) c n = \C E \, <f> = <j>i - <f>2 - 4>3 and, 
si, S2 and S3 are the signs of —Cgw and Ce, respec- 
tively. The coupling coefficients c±, c 2 and C3 are made 
unity by the renormalization 





■+#! 






■+ # 2 




V>3 " 


■+ #3 


= ^CiC 2 1p3 



(83) 
(84) 
(85) 

and we also note that we may, with no loss in generality, 
assume S2 = S3 = 1 (it is only the sign of S1S2 = S1S3 
that matters because an overall sign can be removed by a 
renormalization) . The general solution can then be found 
in litterature, e.g. in Ref. [l£j. The strongest coupling 
occurs for cos0 = ±1 and in this case there are some 
particularly simple solutions that are listed below (recall 
that #2 = #3). 

a) EMW to GW conversion with positive EMW energy 
si = -1, #i(0) = 0, 2,3(0) = # and cos0 = -1 

#1 = #tanh(#t) #2,3 = #sech(#t) (86) 

b) EMW-GW instability with no initial GW 

si = 1, #i(0) = 0, #2,3(0) = # and cos0 = 1 

#1 = #tan(#<) #2,3 = #sec(#£) (87) 

c) EMW-GW instability with equal initial amplitudes 

s x = 1, # x (0) = #2.3(0) = # and cos0 = 1 



# 1 = (#- 1 -i) 1 # 2 ,3 = (# _1 -0 1 
d) EMW-GW interaction with decaying amplitudes 

si = 1, #i(0) = #2,3(0) = # and cos0 = -1 



#1 = #~ x +t) 



#2,3 = (# _1 +f) 1 (89) 
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Also conversion from GW to EMW energy occurs (given 
that there is a small initial EMW perturbation) the solu- 
tion is somewhat more complicated than the ones listed 
above [uj- Note that ^i(O) = * and ^2,3(0) = only 
has the trivial solution ^1 = W and ^2,3 = 0. The so- 
lution a) corresponds to an EMW being converted into a 
GW on the characteristic time scale 7r/2\E r . The solutions 
b) and c) are explosive solutions, i.e. the wave ampli- 
tudes diverges on a finite time t — and t = 1/^, 
respectively. These explosive instabilities can only occur 
for si = s 2 = S3, which is equivalent to saying that the 
GW and EMW energy must be of different sign. The ex- 
plosive instability thus relies on the existence of negative 
energy EMWs and it is the free energy connected with 
the background drift flow that feeds the instability. We 
return to these interesting cases in section lYlII The so- 
lution d) also involves a negative energy EMW, but in 
this case both the EMW and the GW amplitudes tends 
to zero as (VP -1 + 1) -1 . This means that wave energy is 
converted into free energy of the background state, i.e. 
acceleration of the drift flow. As for the cases when the 
wave phases are not fulfilling cos0 — ±1 the behaviour 
is similar 25] but with somewhat larger characteristic 
times. 

We now explore the possibility, indicated by the so- 
lution d), that the waves may interact to accelerate the 
background drift flow. For the wave perturbations we 
are considering here, the fluid equation of motion 117|) 
has the property Vidtjv 1 = g z (where i = 1,2,3). This 
implies that £ = 7(1 — v z ) is an exactly conserved quan- 
tity, as is seen by making the computation 

dti = da - d t jv z = v^av 1 -g z = (90) 

Making use of this, the z-component of the fluid momen- 
tum equation (Eq. I|17|)l may be rewritten as 

t v z = ?f£,A 2 (E 2 h* - E* 2 h) 

= -oj£A 2 \E\ 2 \h\cos((t)) (91) 

where the phase angle <p has the same meaning as in Eqs. 
(|80|I - H82|I . From this equation it is clear that the drift flow 
may be accelerated or decellerated in the ^-direction by 
the interacting waves. Maximum acceleration occurs for 
cos = — 1. In the pump-wave approximation, where 
the waves are assumed energetic and can be taken to be 
of constant amplitude (an approximation breaking down 
when the energy of the drift flow becomes comparable 
with the wave energies), the drift flow v z grows linearly 
in time as (for cos</> = — 1) 

v z =v z (0)+uj£A 2 \E\ 2 \h\t (92) 

and we point out that for relativistic EMW amplitudes 
A 2 \E\ 2 = \v ± \ 2 m 1. 

VII. SUMMARY AND DISCUSSION 

In this paper we have studied the interaction between 
GWs and EMWs in magnetized plasmas for perpendicular 



and parallel propagtion with respect to the background 
magnetic field. The wave evolution equations were de- 
rived in section IIVI and [V] An important point here 
which, as far as we know, has not been fully recognized 
previously is that, given the high frequency approxima- 
tion, the GWs can be taken to be in the TT gauge even in 
the presence of matter. Consequently, the deduction of 
the GW source term associated with the EMW perturbed 
energy-momentum tensor is simplified, see Eqs. (|25|l and 
(|26|l . The GW, on the other hand, produces effective 
currents in the Maxwell equations and an effective grav- 
itational force on the plasma, resulting in a source term 
in the EMW evolution equations and l|48fl. Further- 
more, it was shown in section Ivl that, provided there is 
a velocity drift in the background state, the wave energy 
may be negative [10( for EMWs with phase velocity equal 
to c = 1 (when propagating parallel to the background 
magnetic field). 

Resonant interaction (involving a single EMW) may 
occur when the frequencies matches and the phase veloc- 
ity of the EMW coincide with that of the GW, which to 
a very good precision equals c = 1 in the high-frequency 
approximation. The wave interaction equations govern- 
ing the resonant interaction process were derived and an- 
alyzed in section IvTl for the cases of perpendicular (high- 
frequency extraordinary EMWs) and parallel propagation 
(finite amplitude circular polarized EMWs) with respect 
to the background magnetic field. In both cases the in- 
teraction equations were shown to imply conservation of 
total wave energy. 

i) Perpendicular propagation. The wave interaction 
equations for this case are Eqs. I|67[l and (|68[) . The so- 
lution, given by Eqs. (|69fl and J7DJ, reveals that the 
conversion rate |26| is of order tp = (K/2) 1 / 2 Lu c oj^ 1 B( y 
The effect of the plasma is to diminish the coupling 
strength. The conversion rate becomes ip = ( K /2) 1//2 -B(o) 
for strong magnetic fields, i.e. u> c 3> w p (similar to hav- 
ing no plasma present), whereas for weak magnetic fields, 
lo c <C u) p , the conversion rate is a factor w c /w p smaller. 
Although ip is typically small, GWs may still generate 
EMWs with significant amplitude. For small ipt, the 
EMW are converted into GWs as B = ^u! p u; 2 u;^ 2 B^hot. 

ii) Parallel propagation with positive EMW energy. 
The solutions of the wave interaction equations 172|) and 
(|73|l were found by reformulating them as a three-wave 
interacting system (with two identical EMWs) for which 
the general solutions are well known. From the solution 
(|86|l it follows that the conversion rate [26| for conversion 
from EMW to GW is * = ^C G wC E \E ini \, where E ini is 
the initial electric field amplitude. The conversion rate 
increases with increasing initial amplitude and with in- 
creasing drift velocity (through 7). In order to have more 
transparent formulas we consider a special case, namely 
low- frequency waves, such that u> <C u> C i,u) ce . The con- 
version rate can be formulated as 

*= {CGwClC-^E^mY (93) 
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and in the case of low-frequency waves 

C gw = 
C E = CV- 1 
E E =V\E ± \ 2 

where 

2 2 
C = 7T> -4(1 -U 2 2 , P = l± — > 

and the final GW amplitude is given by hfi na i — 

Hi) Parallel propagation with negative EMW energy. 
As was remarked previously, the EMW energy density 
(given by Eq. I)63|)) may be negative. Such waves are 
perturbations with the property that the energy of the 
perturbed system is less than for the unperturbed sys- 
tem [23|. As is clear from FIG 1, this occurs most com- 
monly for low- frequency waves in electron-ion plasmas, 
i.e. uj <C "feUJpe- The wave energy conservation then 
allows simultaneous growth of the GW and EMW ampli- 
tudes. The physical interpretation is that the background 
flow is unstable with respect to this type of perturbations 
and the energy associated with the flow is converted into 
wave energy. Given that there is an ion flow in the di- 
rection of wave propagation and a large enough plasma 
frequency, this instability will always produce waves with 
a unique EMW frequency determined by Eq. and 
GW frequency given by the matching condition (|66[) . It 
should be noted that this is a nonlinear instability and 
would not have been found using conventional linear sta- 
bility analysis. The solutions l|87|l and (|88|) shows that 
the background is explosively unstable with respect to 
GW-EMW perturbations, the amplitudes reaches infinity 
in a finite time t ~ "J -1 (the interaction equations are of 
course invalidated before this time is reached). The solu- 
tion (|89l) is an example of the inverse process, where both 
waves decreases in amplitude. This means that wave en- 
ergy is converted into free energy of the background state, 
i.e. the drift flow is accelerated. An evolution equation 
for the background flow was also derived, see Eq. (|91|l . 
The response of the plasma on the interacting waves is an 
acceleration (if cos(f> > 0) or deceleration (if cos</> < 0). 
The conclusion is that GWs and EMWs can interact to 
produce a drift flow along the background magnetic field 
with initial linear growth rate as large as ury(l — w z )|/i|. 
The magnitude of the produced drift flow can be esti- 
mated from the initial wave energies and the time scale 
for all the wave energy to be converted into drift flow is 
at least of the order given by Eq. i|89|) . t ~ , E r_1 . 

Finally we discuss the relevance of the results presented 
in this paper to various astrophysical and cosmological 
processes. The phenomena studied in this paper involve 
energy transfer between GWs, EMWs and a background 
flow. This could in principle be an important feature 
in any scenario involving magnetized plasmas, energetic 
GWs and EMWs and/or strong background plasma flows, 



e.g. supernova explosions, gamma-ray bursts, jets of ac- 
creting condensed objects, phenomena in the vicinity of 
neutron stars and in the primordial fluctuations of the 
cosmological plasma. The energy conversion must take 
place within reasonable times and/or volumes for the 
considered processes to be of significance. As a spe- 
cific numerical example let us first consider graviational 
to electromagnetic conversion for perpendicular propaga- 
tion. For the case of GWs from a binary pulsar close 
to collapse (corresponding to a wavelength of, let us 
say, A ~ 10 5 m) we may have h ~ 0.001 a few wave- 
lengths away from the pulsar. A rough estimation based 
on Eq. H69|) then gives a distance L ~ X/h ~ 10 8 m 
for the induced magnetic field to be comparable to the 
unperturbed field |23|. A significant transfer of GW- 
energy is not realistic in this example, however, since 
the background energy density is too low. As another 
example, let us consider EMW to GW conversion for 
parallell propagation and positive wave energy in the 
low-frequency regime. Let us take Bq ~ 10~ 3 T and 
no ~ 10 18 m~ 3 . Naturally, the beam velocities are as- 
sumed to fulfill H58fl . and the wave frequency is assumed 
to obey w«w cl ~2x 10 5 s _1 . Using Eq. (|9"5|) the char- 
acteristic time scale for energy conversion then becomes 
t ~ v]/ -1 ^ io 7 (SE^ini)^ 1 ^ 2 s : where the electromagnetic 
wave energy density should be given in Si-units. Appar- 
ently we need wave energy densities £E,ini ~ 10 14 J/m 3 
to get conversion within 1 second. Note that EMW en- 
ergy densities in the laboratory can be magnitudes larger 
than this [29j and the value is also modest as compared to 
some astrophysical events, like gamma-ray bursts. How- 
ever, we note that the required wave energy content is 
still very large, since the EMW must be distibuted in a 
volume not less than a cube with a side of 1 light second, 
for efficient GW conversion to take place. 

The model we have presented is idealized in several 
ways, e.g. we have assumed a uniform background and 
monochromatic waves. The model can be extended to 
include weakly nonuniform backgrounds, where the in- 
teraction can be treated using mode conversion theory, 
see e.g. Ref. [27j and Refs. therein. It should also be 
noted that the background drift flow is in general also un- 
stable with respect to other types of perturbations, e.g. 
ion-accoustic perturbations which typically may have a 
higher growth rate. With a more realistic velocity distri- 
bution among the particles, including a thermal spread, 
the conditions and growth rates for various instabilities 
change. It is an open question whether GW-EMW insta- 
bilities (of the type presented in this paper) can be the 
only ones (or the dominating ones) in certain regions of 
parameter space. We have also treated the gamma fac- 
tors and the background flow as constants in time. This 
is not a valid approximation for the long term evolution 
when the change in energy becomes comparable to the 
initial energy of the source. An exception to this case is 
interaction between waves with "moderate" positive en- 
ergies in the presence of a highly energetic background 
flow, such that (1 — v^)~z and v z are approximately 



12 



constant for all times. Variations in 7 ad v z can, as these 
occurs in the dispersion relation Q57f) . also affect the fre- 
quency of the EMW, thereby invalidate the assumed per- 
fect resonance and thus reducing the coupling strength. 

In conclusion, energy transfer between GWs, EMWs 
and a background plasma flow should occur also in less 



idealized situations than considered in this paper, but 
in a more involved way. An improved model that cor- 
rectly describes the long term evolution of the interac- 
tion between a GW, an EMW and the background flow 
in an inhomogeneous background is a project for future 
research. 
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